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00 , We develop a new interpolation scheme, based on harmonic inpainting, for reconstructing the 

' cosmic microwave background CMB temperature data within the Galaxy mask from the data outside 

the mask. We find that, for scale-invariant isotropic random Gaussian fluctuations, the developed 
algorithm reduces the errors in the reconstructed map for the odd-parity modes significantly for 
P3 ■ azimuthally symmetric masks with constant galactic latitudes. For a more realistic Galaxy mask, 

^ ' we find a modest improvement in the even parity modes as well. 
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PACS numbers: 98.80.-k, 98.70.Vc, 04.25. Nx 



INTRODUCTION 



After the first release of the cosmic niicrovirave background (CMB) anisotropy data from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) in 2003, various types of "anomalies" in the CMB temperature anisotropy 
on large angular scales have been reported [E 0, 1 i, S i, S B i, [H] • The origin of these anomalies - whether 
' they are cosmological, statistical fluke, or something else - is unknown pH [T^. Tisl. [Tsl. IT^ [TtI ITU . 
I— ' . Recent analyses based on either a Bayesian or frequentist approach show that the statistical significance level 
^ of some of the large- angle anomalies is sensitive to the treatment of the Galactic sky cut^, i20|. The lack of 
, robustness mainly comes from the fact that we cannot use the full sky information, but can only use a part of 
• the sky outside the Galaxy mask. There is a way to construct a full sky map that minimizes the foreground 
I contamination, e.g., the Internal Linear Combination (ILC) technique pll. [2^: however, we cannot rule out a 
lO ■ potential residual foreground contamination in the ILC map, especially in the region that is very close to the 
' Galactic plane, due to a limited number of sky maps at different frequencies. 

An alternative approach is to reconstruct the temperature data within the Galaxy mask from the information 
' available outside the mask. So far, various types of methods such as the direct inversion, the Wiener filtering, 
OO , and the power equalization filtering have been proposed for reconstructing the CMB anisotropy on the cut 
' sky [2^. The Gibbs sampling provides a way to compute an ensemble of the Wiener filtered maps with a 
statistical sample of power spectra fit to a given input sky map, which can be used to estimate the best- fitting 
Wiener filtered map and the corresponding uncertainties; see [25| and references therein. All these methods are 
based upon a linear transformation of the expansion coefficients of spherical harmonics on the cut sky, a'f^, 



into the full sky coefficients, o,f^^\ which are suitable for nearly all-sky data with a small sky cut. However, it 
has been shown that the methods do not work properly for maps with a large sky cut|26i, i27i|. 

In order to reconstruct the data within the Galaxy mask, one needs to assume a prior on the properties of 
the data. For instance, one may require the temperature anisotropics and their absolute values of the gradient 
to have a Gaussian distribution. Then, one can find an optimal solution for reconstructing the data within the 
masked region from available information outside the mask. Such an operation is called an "inpainting" , which 
has been used by skilled museum or art workers for restoring damaged photographs, films, and paintings. In 
recent years, various automatic inpainting algorithms based on partial differential equations or the variational 
principle have been proposed and used for automatic image restoration and removal of occlusions [28l [29l. Isoj . 

In this paper, we formulate an algorithm based on harmonic inpainting for reconstructing smooth CMB data 
within the Galaxy mask. In Sec. II, we develop a numerical scheme for implementing harmonic inpainting on 
a unit sphere based on the boundary element method. In Sec. Ill, we use Monte Carlo simulations to estimate 
the errors of the inpainted signal for azimuthally symmetric sky cuts as well as for a realistic non-symmetric 
cut. In Sec. IV, we summarize our results. 
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II. HARMONIC INPAINTING 



Let X ~ x{9,(j)) be a unit pointing vector on the sky towards a given direction, {0,(j)), and mo(x) be an 
observed temperature fluctuation outside the Galaxy mask, x e Z), where D denotes the region within the mask 
(an open connected set) and D outside the mask (the comphment of D). We assume that the boundary of D, 
dD, is smooth. We wish to reconstruct a smooth temperature fluctuation on the fufl sky, or a best inpainting 
w(x), from the data outside the mask, Uo(x). 

We reconstruct the full sky map, w(x), by locally minimizing the following "energy", 

n[u]= f [X{u~uo)^ + \Vu\^]^dv^, (2.1) 

where A is a positive constant, gij is the Riemannian metric tensor of a sphere, and dv is the infinitesimal 
Euclidean volume. The first term and the second term in the r.h.s of Eq. (|2.ip represent faithfulness to 
the imperfect data and the regularization penalty, respectively. The constant. A, controls a tradeoff between 
faithfulness to the imperfect data and smoothness of the fluctuations. This procedure is optimal when the 
absolute values of the gradient, |Vu|, and the difference between u and uq are Gaussian distributed^. The 
Euler-Lagrange equation of (|2.ip is given by 

- Au + A(-u - uo) = 0; A>0 for x ^ D, A = for xeD. (2.2) 



In order to solve Eq. p.2p . we also need a solution for the derivative of u with respect to the unit vector n 
normal to the boundary of the mask. 



du / du 
dn \ dx^ 



(2.3) 



where x G dD for which u and q are continuous. 

The minimizer of Ea. (j2.ip is called the harmonic inpainting denoizing of mq. Expanding a temperature 
fluctuation in real spherical harmonics as u = 'YliO-imYim and uq — 'YlibimXim, the solution of Eq. (j2.2p for the 
region outside the mask, x € l), is given by 

«?r(A) = ^(7^^.™, (2.4) 

where Q°^*'^ ^'"^ ^^^^ expansion coefficients for the solution outside the mask. Eq. (|2.4p implies a suppression of 
modes with I ^ ^/X. In the limit of A ^ 1 ^, the minimizer of Eq. (|2.2p is called the real harmonic inpainting. 
The fiuctuation, m(x), satisfies the following boundary integral equation (see [3l| for further detail). 



u(x) 



/ G(x, y)g(y) Vffrfy - / ^(x, y)u{y)^dy = 0, x e dD, 

JdD JdD 



(2.5) 



where G is the Green's function of Laplacian A, and H is the normal derivative oi G, H = dG/dn. By dis- 
cretizing the boundary, dD, into 2N elements, Tj, and approximating u on dD by some low-order polynomials, 
Eq-daSl) yields 

[H]{u}^[G]{q}, q^^, (2.6) 

where {u} and {q} are iV-dimensional vectors consisting of the boundary values of Ui and their normal deriva- 
tives, qi, and [H] and [G] are x A^- dimensional coefficient matrices which are obtained from the integration 
of the fundamental solution, G(x, y), and its normal derivatives, iJ(x, y). For constant elements 

M(xj) — — Const. , 9(xj) = q^ = Const, on Tj, (2-7) 



^ This is a reasonable requirement for Gaussian random fluctuations as each component of the gradient is also a Gaussian random 
variable. 

^ This corresponds to the case where the Gaussian noise of the fluctuation u — uq outside the mask is negligible in comparison 
with the expected absolute values of gradients |Vu|. For instance, from the signal-to- noise ratio, one can estmate that A > 10'' 
for the WMAP data at I < 15. 
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the elements of the matrices are 



where 



Hii — 



Hij i ^ j 



Htj = I |j^(x,,yj) y/gdy, Gy = / G(x,, y^) ^gdy. (2.8) 

In terms of the normal derivatives {q\ that are derived from Eq. (j2.6p . the reconstructed function u within 
the masked region, D, is given by 

u{^) = ~ [ G(x,y)<z(y)Vgrfy+ / H{x,y)u{y)^dy, xeD. (2.9) 

JdD JdD 

III. ERROR ESTIMATION 

A. Application to azimuthally symmetric masks 

In order to check the robustness of the algorithm developed in section II, we carry out Monte-Carlo simulations 
of random Gaussian fluctuations on a unit sphere, Ug. We generate Monte Carlo realizations from a scale 
invariant power spectrum (n = 1), 

Ci - {\airn\') oc (3.1) 

which describes the large-angle CMB temperature fluctuations in the Einstein-de Sitter (EdS) universe. 

We first consider azimuthally symmetric masks with constant galactic latitudes, i.e., we mask region below 
given galactic latitudes, |6| < 9q. For isotropic random fluctuations, modes with I > tt/Oq are not affected by 
the azimuthal cut very much (except for the overall normalization). Therefore, in the following we shall use 
only multipoles with I < Ic, where Ic — tt/Oq. From Eq. (|2.4p . one can see that this choice of cutoff in harmonic 
space is roughly equivalent to choosing the control parameter. A, such that A = (tt/^q)^ = 'c- The number of 
elements on the boundary necessary for numerically solving the integral equation (j2.5|) depends on the size of 
the mask. For instance, we discretize the boundary of the mask into 80 linear elements for 9q = 20°. 

To measure the difference between the reconstructed fluctuation, u(x), and the original one, %(x), we use 
the L^-norm defined as 

{a[u\f^( {u[e,^]fd^l^Y.^afJ,)\ (3.2) 

where (Kl is the infinitesimal area element for a surface of unit sphere, and a|^^j is the spherical harmonic 
coefficients of the reconstructed map. To reduce the sample variance, we generate 1000 realizations, each of 
which is masked by the azimuthally symmetric mask and is reconstructed by the developed algorithm. We then 
calculate the mean relative errors defined as 

/Aax2_ / /g,(ug[g,0]-u[g,^])^dr> \ 

\ a) - \ /g,(u<,[0,^])2df] /' ^ • ^ 

where Mg(x) and u(x) are the underlying Gaussian fluctuation and the inpainted (reconstructed) fluctuation, 
respectively. For comparison, we also compute the relative errors for the fluctuations with the "naive" ansatz 
in which w = within the masked region and u = uq outside the mask. 

As one can see in Fig[l] some of the missing features within the masked region are reconstructed by the 
inpainting. Improvements are more conspicuous for odd multipoles. 

In Table HI we compare the mean relative errors, (Aa/a)'^, at Z = 1, 2, 3, 4, and 5, for two azimuthally 
symmetric masks, \h\ < 20° and 30°. Compared with the naive cases, our developed algorithm based on harmonic 
inpainting significantly improves reconstruction of odd multipoles. However, for azimuthally symmetric masks, 
we find that the inpainted even multipoles are identical with those in the naive case. 




FIG. 1: Contour maps of one Gaussian realization (top), those with an azimuthally symmetric mask \b\ < 20° (middle), 
and the inpainted maps (bottom) in MoUweide projections. The left figures show co-added multipoles (from / = 1 to 
I = 5) whereas the right figures show co-added odd multipoles (/ — 1,3,5). The boundaries of the mask are shown in 
white lines. 



\b 


< 20° 


I 


1 


2 


3 


4 


5 


naive 


0.23 


0.21 


0.55 


0.39 
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inpainted 


0.13 


0.21 


0.31 


0.39 


0.32 


\b 


< 30° 


I 


1 


2 


3 


4 


5 


naive 


0.39 


0.37 


0.76 


0.58 


0.62 


inpainted 


0.27 


0.37 


0.50 


0.58 


0.48 



TABLE L Mean relative errors, (Aa/a)^, of multipoles for azimuthally symmetric masks with constant galactic latitudes 



This can be explained as follows. The Green's function on a unit sphere is given by G(x, y) = —Qq{x ■ y)/27r, 
where 2(5o(x • y) = ln[(l -f x • y)/(l — x • y)]; thus, we have G'(x, y) = G'(— x, — y). In other words, the Green's 
function is invariant under the spatial inversion. Therefore, if the boundary of the masked region is invariant 
under the spatial inversion, the order of the 2A^-dimensional element matrix [G] is reduced to iV, which kills 
the degree of freedom for even-parity modes. 

The sky cut biases the estimation of the angular power spectrum when it is computed naively as Ci(est) = 

X]m('^]m)^/(2^ + 1)' which is often called the "pseudo C" (s^]- The effect of mask needs to be deconvolved 
using, e.g., the MASTER method [s^. Here, we do not use the MASTER method, but study how much the 
bias is reduced when we use the naive C; estimation and the harmonic inpainting. In Figure[2]we show the ratio 
of the estimated C; divided by the input Ci, averaged over 1000 Monte Carlo simulations, for three azimuthally 
symmetric masks, |6| < 10°, 20°, and 30°. We find modest improvements at odd multipoles, whereas no 
improvements are found at even multipoles as they are identical to those in the naive case. 
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FIG. 2: The estimated (pseudo) power spectrum, Ci(est) = ^ \au 
\h\ < do = 10°, 20°, and 30°, divided by the input Ci oc {l{l + 1)) 
1000 Monte Carlo reahzations, and plotted for the inpainted (reconstructed) fluctuations (solid lines) and for the naive 
ones (dashed lines). 



P/(2i + 1), for azimuthally symmetric masks with 
^. The estimated Ci(est) have been averaged over 



B. Application to a realistic Galaxy mask 



As we have seen, the developed algorithm based on harmonic inpainting fails to reconstruct even multipoles 
when the Galaxy mask is azimuthally symmetric. However, the shape of the Galactic foreground emission 
region is not perfectly azimuthally symmetric, and thus the corresponding Galaxy mask does not respect parity 
conservation. Therefore, we expect an improvement for even multipoles as well relative to the naive estimations, 
once realistic Galaxy masks are used. 

To confirm this expectation, we have carried out Monte-Carlo simulations with a realistic Galaxy mask. In 
this study we also use a realistic sky signal as well: we have generated 1000 realizations of temperature maps 
from a ACDM model with cosmological parameters given by {n,flm,^A) — (1,0.24,0.76). As for the Galaxy 
mask, we use the WMAP's KpO mask[2T]. 

Since the shape of the KpO mask is fairly irregular at the scale of pixels, the definition of the border is not 
so trivial like the azimuthally symmetric mask in Fig. 1. For instance, some observed pixels close to the mask 
are surrounded by two or more pixels, some others are completely surrounded by masked pixels, making the 
choice of the border arbitrary and unstable. To avoid this we have built a smoothed version of the KpO border. 
Our procedure consists of going (for a fixed sampled with a constant interval) along the coordinate 9 from 
Galactic north to south until we find a masked pixel belonging to the KpO mask. Once we have found this 
pixel, we would go back by a step of 10 arcminutes and choose that pixel living on the border. The final set of 
(&i, 4>i) defines our border. We found that this procedure regularizes the KpO border. 

We choose the cutoff angular scale to be Ic = 15, which corresponds to the mean angular size of the mask 
in the polar direction. The boundary of the KpO mask are discretized into 80 linear elements {6i,(l)i) with a 
constant interval 27r/40 (Fig[2]). 

As in the case of the azimuthally symmetric mask, some of the missing features within the KpO mask are re- 
constructed by the inpainting. The improvements are more conspicuous for odd multipoles than even multipoles 
(FigSl). 

We show the mean relative errors for the naive estimation and those for the inpainted map with the KpO 
mask in Table IL We find that our algorithm improves reconstruction of both the odd and even multipoles, but 
the improvement for the even multipoles is modest (3 to 10%). 

In order to find the polar-angle dependence, we compute the mean relative errors for each (/, m) mode, defined 
by 




where af^ and a;,„ are real expansion coefficients for the inpainted map and the original map, respectively. 

We find that significant improvements come mainly from odd-parity modes with m = ±1 (Fig. [5]). Lower 
multipoles are improved more by inpainting than higher ones. We also find that the biases in the pseudo C; are 
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FIG. 3: Meshes on the boundary of the KpO Galaxy mask, smoothed on scale ~ 10 arcmin. 



/=1-6 /=H3f5 /=2+4+6 




FIG. 4: Contour maps of one Gaussian realization (top), those with a KpO mask (middle) and the inpainted maps 
(bottom) in MoUweide projections with galactic coordinates. The left figures show co-added multipoles (from / = 1 to 
/ = 6) whereas the middle and right figures show co-added odd multipoles {I = 1, 3, 5) and co-added even multipoles 
[l = 2, 4, 6), respectively. 

reduced for both the odd and even multipoles (Fig. [S]). 

IV. SUMMARY 

We have shown that our new algorithm based on harmonic inpainting offers a way to reconstruct the missing 
CMB temperature anisotropy data within the Galaxy mask, provided that the primordial fluctuation is Gaussian 
with a scale- invariant {n — \) spectrum. The method reduces the errors for multipoles with an odd-parity (odd 
I) significantly. While the method fails to reconstruct multipoles with an even parity (even I) for azimuthally 
symmetric masks, it can reduce the errors for multipoles with an even parity for a more realistic, non-azimuthally 
symmetric mask, such as the WMAP's KpO mask. As we have shown, our method is useful for studying the 
dependence of the CMB at low multipoles given by the azimutal parameter m of the harmonic coefficients, 
which is relevant to confirm or confute some of the anomalies at low multipoles found in the WMAP data. 

Our new algorithm does not require any specific form of the power spectrum as a prior for reconstructing 
the missing data as long as the fluctuations obey a Gaussian distribution. Even when fluctuations are not 
Gaussian, our algorithm is expected to work as long as the fluctuations are sufficiently smooth and the gra- 
dients approximately obey a Gaussian distribution. Therefore, our method may be used to reconstruct other 
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naive 


0.096 


0.11 


0.28 


0.23 


0.35 


0.28 


0.35 


0.30 


0.33 


inpainted 


0.053 


0.10 


0.15 


0.21 


0.21 


0.27 


0.24 


0.29 


0.25 



TABLE II: Mean relative errors (Aa/a)^ of multipoles for the KpO Galaxy mask. 




FIG. 5: Mean relative errors (Aa/a)^j for the KpO Galaxy mask. The solid lines and the dashed lines show the 
inpainted and naive cases, respectively. 

observational data that is sufficiently smooth such as radio or infrared maps on a cut sky (a similar technique 
that can incorporate texture as well as the smooth part is studied in [s^). 

Our analysis implies that low multipoles with an even parity are more likely to suffer from the effect of sky 
cut compared with those with an odd parity. In other words, even parity modes on large angular scales cannot 
be precisely reconstructed without assuming specific forms of the angular power. Therefore, the statistical 
significance of the anomalous features involving the CMB quadrupole may be much lower than claimed, as 
recent studies suggest 19, 20]. 

The feature of the inpainted 5 year WMAP data will be explored in the forthcoming paper (Inoue & Cabella 
2008). 
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